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ABSTRACT 

V»> , We simulate the formation and chemodynamical evolution of 124 elliptical galaxies 

■ by using a GRAPE-SPH code that includes various physical processes associated with 

the formation of stellar systems: radiative cooling, star formation, feedback from Type 

t-H ■ II and la supernovae and stellar winds, and chemical enrichment. In our CDM-based 

' scenario, galaxies form through the successive merging of sub-galaxies with various 

<^> masses. Their merging histories vary between a major merger at one extreme, and a 

^^O ■ monolithic collapse of a slow-rotating gas cloud at the other extreme. We examine the 

physical conditions during 151 merging events that occur in our simulation. The basic 
processes driving the evolution of the metallicity gradients are as follows: i) destruction 
by mergers to an extent dependent on the progenitor mass ratio, ii) regeneration 
when strong central star formation is induced at a rate dependent on the gas mass 
of the secondary, iii) slow evolution as star formation is induced in the outer regions 
through late gas accretion. We succeed in reproducing the observed variety of the 
CI.' radial metallicity gradients. The average metallicity gradient A log Z/ A log r ~ —0.3 

q . with dispersion of ±0.2 and no correlation between gradient and galaxy mass are 

?-h ' consistent with observations of Mg2 gradients. The variety of the gradients stems from 

, the difference in the merging histories. Galaxies that form monolithically have steeper 

C$ ■ gradients, while galaxies that undergo major mergers have shallower gradients. Thus 

| merging histories can, in principle, be inferred from the observed metallicity gradients 

• i-h , of present-day galaxies. The observed variation in the metallicity gradients cannot be 

■ explained by either monolithic collapse or by major merger alone. Rather it requires 

a model in which both formation processes arise, such as the present CDM scheme. 

CS | 
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1 INTRODUCTION 

The internal structure of galaxies, spectrophotometric, 
chemical, and dynamical properties of various locations 
within a galaxy, is closely related the processes of galaxy for- 
mation and evolution. Stars in a galaxy are fossils; the star 
formation and chemical enrichment history of the galaxy 
are imprinted on their kinematics and chemical abundances. 
The SAURON pr oject with William Herschel Telescope 
jBacon et alj UoOl) is providing the wide-field mapping of 
the kinematics and stellar population, which will certainly 
give stringent constraints on the galaxy formation and evo- 
lution. Multiobject and Integral Field Spectrographs are be- 
ing developed also on 8-10m ground-based telescopes, which 
can give the time evolution of such internal structure. To 
derive the physical processes from such observational data, 

* E-mail: chiaki@MPA-Garching.MPG.DE 



it is necessary to construct a realistic model, i.e., a three- 
dimensional chemodynamical model, and to compare the 
theoretical predictions with such observational data. 

Theoretical approaches have been perf ormed in many 
ways: 1) The one-zone model (e.g ., iTinslevI ll98(J : 
lArimoto fc YoshiillT987l iMatteuccil 11999) played an impor- 
tant role in constructing the basic evolutionary scenarios 
of galaxies. In this model, the star formation history of a 
galaxy is constructed by carefully comparing the model pre- 
dictions with the mean photometric and chemical properties 
of observed galaxies. However, because of the simplified as- 
sumption that all matter in a galaxy is well-mixed instan- 
taneously, the one-zone model provides no idea of the inter- 
nal structure of galaxies. 2) The s emi-a n alytic mod e l (e.g . , 
iKauffmann. White fc Guiderdonil Il993l : ICole et alJ Il994ll . 
which is based on the Press-Schechter theory, can provide 
the mass function of dark halos, their survival timescales, 
and the merging rates. By adopting empirical laws to de- 
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termine the stellar mass and the mass of heavy elements, 
and by introducing simple rules to determine the galaxy 
morphology (i.e., an elliptical galaxy forms from the ma- 
jor merger of spiral galaxies), the semi-analytic model could 
reproduce some correlations among global properties (e.g., 
the Tully-Fisher relation and color-magnitude relation) and 
some constraints on the number of galaxies (e.g., the lumi- 
nosity function and number counts). However, there seems 
to be some difficulties t o explain the number evolution o f 
elliptical galaxies (e.g., iBenson. Ellis fe Menanteaul l20o3) . 
The observed information on the internal structures remains 
untouched . 3) Numerical simulations of dissipa tionless sys- 
tems fe.g.. lToomre fc Toomrell972l ; lw"hitell978l) provide the 
interpretation of the interaction of galaxies and the inter- 
nal structure of a galaxy. The gas dynamics in three di- 
mensions was included in such numerical simulations (e.g., 
iHern quist fc Kat3ll9 89). and then star formation, feedback 
jKgt zl ll992l) . andchemical enrichment l|Steinmetz fc Miilleil 
11994) were included. It is now possible to predict the spatial 
distributions of gas, stars, and heavy elements in a galaxy. 
However, the comparison with the observations has not been 
fully attempted yet. The reason is that it takes long time to 
calculate the evolution of one galaxy with enough resolution 
to predict such distributions. To reduce the calculation time 
and to improve the resolution, a variety of methods have 
been invented. 

The Smoothed Particle Hydrodynamics method (SPH) 
is widely used to calculate three dime nsional hy- 
drodynamics with a Lagr a ngian scheme llLucvl Il977t 
lOingold fc Monaghanl Il977t iMonaghanl 11993 for a re- 
view), and has been applied to many astrophysical prob- 
lems th at have large density contrasts; formation of 
galaxies jHernnuist fc Ka^td \lwl iNavarro fc White] ITool 
Stcinmetz & Miillcr |l994 |) and a cloud -cloud collision 
IILattanzio et al.l 11985 : lHabe fc Ohtalll992T) . Various codes 
have been developed to combine SPH with collision- 
less particles (i.e., dark matter and star particles; here- 
after an N-body system) using different methods to cal- 
culate gravitational forces; dir ect summatio ns, Particle- 
Parti cle/Particle-Mesh metho ds ()E\3^rd|[l988j) , Tree meth- 
ods jHernauist fc KatJ Il989l : iBenz et alJ ll99Cll) , and the 
method using the special purpose computer GRAPE 
JUmemura et al.lll993l ; lsteinmetzlll996l) . GRAPE (GRAvity 
PipE) is a special purpose computer for efficiently calculat - 
ing gravitational force and potential dSugimoto et 81.1 11990). 
The GRAPE-SPH enables us to simulate the formation and 
evolution of a galaxy with more than 10 4 particles in calcu- 
lation time as short as a few days. It makes it possible to 
simulate many types of galaxies with different initial condi- 
tions, which is crucial to study the formation and evolution 
of galaxies statistically. 

The aim of our study is to put constraints on the 
formation history of elliptical galaxies by comparing the 
observed internal structures of stellar population. To con- 
struct a self-consistent three-dimensional chemodynamical 
model, we have introduced various physical processes as- 
sociated with the formation of stellar systems such as ra- 
diative cooling, star formation, feedback of Type II and 
la supernovae (SNe II and SNe la), and stellar winds 
(SWs), and chemical enrichment. The chemical enrichment 
of SNe la has been recently included in several chemo- 
dynamical models. Among two alternative scenarios of 



the SN la progenitor (e. g.. iKo bavashi et alJ Il9 98l). most 
chemodynamical mod e ls jRaiteri. Villataj£jjavar£o| ]l99d 
ICarraro. Lia fc Chiosil ll998Tl adopted HGTe^rgio^^Remdnil 
Jl983ft 's formulation based on the double-degenerate sce- 
nario. iMosconi et alJ (1200 if) adopted single time delay for 
the SN la contribution, and the parameter corresponds to 
the double-degenerate scenari o. The single-degenerate sce- 
nari o has been introduced in iNakasato fc Nomotol J2003T) 
and iKawa^^^G^bsonl ( 2002 ). It may be useful to note 
that lWooslev fc Weaveil {l995)'s iron yield is too large com- 
pared with the observed abundance ratios in the Milky Way 
Galaxy, and thus their iron yield should be reduced to be 
half. Such modification is always adopted in the one-zone 
models, but is never mentioned in many chemodynamical 
models. We have constructed a realistic model of chemical 
enrichment, excluding the instantaneous recycling approxi- 
mation, including the mass-dependent yields of SNe II and 
the single-degenerate scenario of SNe la. We then solved 
the evolution of slowly-rotating systems that consist of dark 
matter, gas, and stars from various initial conditions to pre- 
dict the spatial distribution of stellar population within a 
galaxy. By comparing the theoretical metallicity gradients 
with the observed ones, we discuss the origin of elliptical 
galaxies. 

How elliptical galaxies form is a long-standing issue as a 
matter of big debate. The regularity in the light distribution 
and the global velocity anisotropy i n elliptical galaxies were 
explained by the violent relaxation jLvnden-Belllll967fr) . Ef- 
fectively dissipationless formation of an elliptical galaxy was 
discussed in various ways; e.g., stars have f ormed prior to the 
beginning of the collapse of the gas cloud (lGotdll973l. 1975) 
or stars have formed slowly in disk galaxie s which subse- 
quently merge to make a spheroidal ga laxy l|Toomrelll977l : 
iMarchant fc Shapirdll977l: lBarneslll988l) . The limit of com- 
pletely dissipationless collapse is amenable to N-body exper- 
iments, which can form the objects that have the observed 
dynamical properties. However, the dissipation during the 
formation is indispensable to explain the photometric and 
chemical properties of elliptical galaxies such as the color- 
magnitude relation, the mass-metallicity relation, and the 
radial metallicity gradients. 

Two competing scenarios of the formation of elliptical 
galaxies have so far been proposed: Elliptical galaxies should 
form monolithically by gravitational collapse of gas cloud 
with considerable energy dissipation (herea fter referred to 
as the monolithic colla pse hypothesis; e.g., iLarsonl Il974bt 
lArimoto fc Yoshiil Il987l) . or alternatively ellipticals should 
form via mergers of gaseous disk galaxies or of many dwarf 
galaxies (here after referred to as the merger hypothesis; e.g., 
| Toomrelll977l: Kauffmann et al. 19 93; iBaugh. Cole fc Frenkl 
Il996t ISteinmetz fc Navarrol 120021) . The merger hypothe- 
sis can be supported by the dynamical disturbances of 
observ ed ellipticals such as shells /ripples and mu ltiple 
" 19901: ISchweizer fc SeitzeiHl992l; see 
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1992), and may easily ex plain the 



morphology-density relati on of galaxies in clusters jDressleJ 
ll980l : lDressler et al.lll99^) . However, elliptical galaxies show 
apparently little evidence for o n-going star formation, the 
bulk of their stars are old (e.g.. [Kodama fc Arim otolll997|: 
[Stanford j^iserihajdtfc Dk^nsonl ^98[ iKo^jamaet^tl] 
Il998t ISilva fc Bothunl Il998ft . The monolithic collapse 
hypothesis assumes that the bulk of stars in ellipticals 
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form during an initial star burst at high redshift, and that 
the star formation is terminated by a supernovae-driven 
galactic wind that expels the left-over interstellar gas 
from galaxies. The galactic wind is supposed to play an 
essential role in injecting heavy e lements into the hot 
intracluster gas l|Ciotti et al.l Il99lf ). and predicts tight 
correlations among global pr operties of galaxies such a s 
the color-magnitude relation feower. Lucev fcjShM^ Jl992j). 
the metallicity-velocity dispersion rel ation jDavies et aD 
1987) , and the fundame ntal plane jDiorgovski fc Davisl 



19871: iDressler et al.lll987l) . (The color-magnitude relation 



could be reproduced also u nder the merger hypothes is 
jKauffmann fc Chariot] Il99gl. but see ICole et al.l |2000|U 



Recent observations of clusters at high red shifts reveal that 
these relationships exist even at z ~ ljDickinsorj 
ISchade. Barrientos fc Lopez-Crudll997l : iKelson et al " 



1996; 



1997 



Stanford et al. 1998), which indicates that the bulk of 
stars in cluster ellipticals forms at the redshift zt > 2.5 — 4 
l|Kodama et alJll998h . 

However, for cluster early-type galaxies, it has been 
argued that the "progenitor bias" is significant; the pro- 
genitors of the youngest low-redshift e arly- type galaxies 
drop out of the sample at h igh redshift <IKauffmanrJll996t 
Ivan Dokkum fc Franxl Il996l. 2001). Thus the evolution of 
field early-type galaxi es is now paid attention as the ob- 
servational constraints. lFranceschini et alJ ]l998l) found that 
HDF-N early-type galaxies are relatively young with the for- 
mation epochs spanning 1 < z < 4, but no evolution of 
field e lliptical galaxies found at l east z < 1 iSchade et alJ 
[ 19991: iBrinchmann fc Ellis! l2000t llm et alJ l2002t s ee also 
baddi. Cimatti fc RenzinilboOOl) . iDrorv et al] J200 ll) found 
mass evolution at 0.4 < z < 1.2. At high redshift, the num- 
ber of field red ellipticals is s maller than expected by the 
monolithic collapse hypot hesis (IZepil997l : lMenanteau et al] 
119991: iBarear et aflll999h . and there may be some global 
evolution that is consistent with the hierarchical cl ustering 
scenario jFontana et al]ll99St l Dicki nson et al 1l2003t but see 
ICimatti et al]l2002l) . ~ 

The radial metallicity gradient gives one of the most 
stringent constraints on the galaxy formation. Numerical 
simulations of the collapse of galaxies including star for- 
mation definitely pre dict strong rad ial gr adients in chem - 
ical enrichment (e.g.. lLarsonlll974al 1975: ICarlberdll984l) . 
whereas the dissipationless co llapse mod els predict no gra- 
dient in chemical enrichment <GotJl97l 1975). During the 
collapse, gas is chemically enriched, flows inward, and forms 
new stars, which form the radial metallicity gradients. The 
metallicity gradients are observed as radial gradients of col- 
ors and spectral line indices (e .g., |Faber||l97j; |Faber| W77l 
iDavies. Sadler fc Peletie3ll993t: iKobavashi fc Arimotolll999] 



for a review). A typical observed metallicity gradient of el- 
liptical galaxies is A log Zj A log r ~ —0.3, which is less steep 
than those predict ed by numeric al simul ations of dissipative 
co llapse (—0.35 in iLarsonl Il974at -1.0 in lLarsonlll975l : -0.5 
m ICarlberdll984Tl . Furthermore, if elliptical galaxies form 
monolithically from a massive gas cloud, metallicity gradi- 
ent should correlate with global properties of galaxies in 
the sense that m ore massive galaxies have steeper gradients 
<Carlberdll984h . 

The observational feature of metallicity gradients is 
complicated and confusing because of a lack of suit- 



able sample of uniform quality. It was shown that 
elliptical galaxies with larger values of the central 
Mg2 (~ 5100A) index tend to have steeper Mgg 
gradients iGorgas. Efstathiou fc Aragon-Salamancal Il990l: 



ICarollo. Xi&nz'^^^^usa^QQ^I^a^aI^^i^OTsSt\^9&i . 

However, Davies et al. (1993) did not find any significant 
correlation between the Mg2 gradient and oo in the sam- 
ple of 13 galaxies. Kobayashi & Arimoto (1999) re-studied 
line-strength gradients of 80 elliptical galaxies by using the 
indices of Mg 2 , Mg b (5177A), Fei (5270A), Fe 2 (5335A) and 
H/3 (4861 A), and found that the metallicity gradients do not 
correlate with any physical properties of galaxies, includ- 
ing central and mean metallicities, central velocity disper- 
sions, absolute B-magnitudes, absolute effective radii, and 
dynamical masses of galaxies. Elliptical galaxies have differ- 
ent metallicity gradients, even if they have nearly identical 
properties such as masses, luminosities, and metallicities. 

This discrepancy could be solved if mergers flatten the 
original gradient. Indeed numerical simulations showed that 
the gradient in a disk galaxy should be halv ed after thre e 
successive mergers of galaxies with similar size llWhitell98Cr> . 
However, according to the dissipationless N-body experi- 
ment, the initial state is not fully wiped out during the vio- 
lent relaxation phase, and N-body particles that were in the 
outer region of a progenitor gala xy are found in th e simi- 
lar location after merging events Jvan Albadalll982h . From 
this point, it has been mentione d that metalli city gradients 
is not reduced by a merger (e.g.. lBarneslll996fi . Simulations 
of both dissipative collapse and mergers leave room for im- 
provements, because essential physical processes such as star 
formation, feedback of supernovae, and metal enrichment 
were not taken into account. 

Here we simulate the chemodynamical evolution of el- 
liptical galaxies based on the CDM picture. In the CDM 
cosmology, the amplitude of primordial fluctuation decreases 
with increasing wavelength, and the formation of structure is 
driven by the hierarchical clustering. Galaxies should form 
through the successive mergings of sub-galaxies with vari- 
ous masses. Contrary to the semi-analytic models, we ex- 
clude the assumption that elliptical galaxies form only from 
the major merger of disk galaxies. Instead we allow various 
merging histories for elliptical galaxies. In some cases, an el- 
liptical galaxy forms by an assembly of gas rich small galax- 
ies, which looks like a monolithic collapse. In other cases, the 
evolved galaxies with little gas merge to form an elliptical 
galaxy. This scenario is the midway of monolithic collapse 
and major merger of disk galaxies. 

In this paper, we first describe the GRAPE-SPH code 
and the modeling of physical processes (J2J- We classify sim- 
ulated galaxies according to their merging histories (" H3.H 
and derive the present metallicity gradients ( H3.2H . In iJ3.3l 
we show that the scatter of the metallicity gradients comes 
from the difference in merging histories, and discuss the ori- 
gin of elliptical galaxies by comparing with the observation. 
In H3.4I we examine the evolution of metallicity gradients 
via merging events that occur in our simulation, and man- 
ifest the dependences on mass ratios of merging galaxies, 
gas fractions, and induced star formation. In |3] we mention 
some future works and possible problems. Our conclusions 
are given in Sj^| 
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2 CHEMODYNAMICAL MODEL 
2.1 Hydrodynamics 

Here we describe the GRAPE-SPH code, which is originally 
written by N. Nakasato (2000), and is highly adaptive in 
space and time by means of individual smoothing lengths 
and individual timesteps. The SPH formulation use d in the 
code is almost the same as iNavarro fc White! <1993h . 

Instead of the continuity equation, the density p of each 
particle is determined from 

(p(r i ))=p i = Y,m j W(r,h), (1) 

where r is defined as r = \n — rj\, and h denotes the sym- 
metrized smoothing length defined as h = hij — i(hi + hj). 
The momentum equation is represented as 

= - E m > ( 3 + 3 + ni >) v ^ r - ^ - ( 2 ) 

The energy equation is represented as 

^ = _^ m /^+l niJ V Bi _ rj ).V 4 W(r,A)+2^,(3) 

where thermal conduction is neglected. We also use the equa- 
tion of state for an ideal gas with 7 = 5/3; 

Pi = (7 - l)piUi. (4) 

The smoothing length h varies spatially, evolves 
with time, and is computed for each particle in every 
times tep. For the kernel W, a spherically symm etric spline 
(Monagha n fc Lattanzid ll985: Mo naghanlll992l) is adopted. 
For the derivative of the ker nel VW, the revised form by 
iThomas fc Couchmanl (^992) is ado pted. For the artificial 
viscou s term Hij, we use the modified iMonaghan fc Gingoldl 
(1983) tensor and the shear free viscosity formulation with 
a = 1.0 and j3 = 2.0 IIBalsaral 119951 : INavarro fc Steinmeta 
119971) . The heating TL and cooling A rates are described in 
the next session. 

In the SPH code, a gas particle interacts with dark mat- 
ter and star particles only by gravity. The gravitational po- 
tential $ is given as 

*(r,) = -G£ (5) 
y/\ri - TjY + e 

where G is the gravitational constant, e is the gravitational 
softening length, which we set e = 0.5 and 1.0 kpc for 
the high- and low- resolution, respectively. The gravity be- 
tween these particles is calculated in direct summation using 
the GRAPE5 MUV (Mitaka Underground Vineyard) sys- 
tem in the National Astronomical Observatory of Japan and 
GRAPE6 of the University of Tokyo. 

We use an individual timestep scheme with the follow- 
ing two steps. First, we compute the timestep of each parti- 
cle At c from the d ynamical criteria using the velocity v and 
the acceleration a jKatz. Weinberg fc Hernauisll996t) : At e 

2 

is the minimum of 3-f and n, where ri is a numerical 

\v\ 'V M ' ' 

parameter, and is set to be rj — 0.5. Then, we derive the 

actual timestep Ai a from the greatest power of 2 subdivi- 
sion of the system timest ep At BYB , which is smalle r than At c ; 
At a = At sys /2 n < At c iNavarro fc Whitelll993ft . The sys- 
tem timestep is a fundamental timestep used to synchronize 



all particles, and is set to be At sys = 2 Myr. The time inte- 
gration of the equation of motion is done using a leap-frog 
method modified for the individual timestep scheme. 

2.2 Physical Processes 

To simulate the formation and evolution of stellar systems 
from gas, we introduce various physical processes into the 
GRAPE-SPH code; radiative cooling, star formation, feed- 
back of SNe II, SNe la, and SWs, and chemical enrichment 
including the mass dependence of SNe II. Here we describe 
the formulation and the assumptions of each physical pro- 
cesses. 

2.2.1 Radiative Cooling 

Radiative cooling is modeled using an equilibrium cooling 
function. If gas is primordial with no heavy elements ([Fe/H] 
< —5), we compute the cooling rates using two-body pro- 
cesses of H and He, and free-free emission (Katz et al. 1996). 
These processes are collisional excitation of neutral hydro- 
gen (H°) and singly ionized helium (He + ), collisional ion- 
ization of H°, He , and He + , standard recombination of 
H + , He + , and He ++ , dielectric recombination of He + , and 
free-free emission. We compute a lookup table, which lists 
(Ti. — h)/n\ as a function of temperature and density. We 
then evaluate net cooling rates at intermediate values with 
cubic spline interpolation. 

For metal-enriched gas ([Fe/H] ^ —5), we use a 
met allicity- dependent cooling function computed with the 
MAPPINGS III software by R.S. Sutherland (MAPPINGS 
HI is the updates of MAP PINGS II that is described in 
ISutherland fc Dopitalll993l) . The included processes are col- 
lisional line radiation, free-free and two-photon continuum, 
recombination, photoionization heating, collisional ioniza- 
tion, and Compton heating. Heavy elements can significantly 
enhance the cooling rate. At T > 10 7 K, Fe group line 
emission processes largely determine the cooling function, 
whereas at lower temperatures the lighter atoms such as 
C, O, and Ne dominate. Although the metallicity effect is 
small at T ~ 10 4 K, the cooling rate with [Fe/H] = is ~ 100 
times larger than that for the primordial gas around T ~ 10 5 
K. The radiative cooling, and hence the star formation rate, 
strongly depends on the metallicity. In the table of cooling 
function that we use, cooling rates are given as functions 
of [Fe/H], and the elemental abundance ratios are set to be 
constant for given [Fe/H] according to the relations in the 
solar neighborhood; Galactic halo stars (i.e., [0/Fe]= 0.5) 
for [Fe/H] ^ — 1, and the solar values for [Fe/H] ^ 0, and 
interpolated the halo and solar values for —1 < [Fe/H] < 0. 
Since the elemental abundance ratios depend on the star 
formation history and the inhomogeneous mixing of SNe II 
and la, the cooling tables depending on [O/Fe] are required 
to increase the accuracy. 

2.2.2 Star Formation 

The t reatment of star formation is similar to that in iKatzl 
(1992), which is widely used in GRAPE-SPH simulations. 
If a gas particle satisfies the following star formation crite- 
ria, a fractional part of the mass of the gas particle turns 
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to a star particle. Since an individual star particle has the 
mass of 1O 5_7 M0, it dose not represent a single star, but an 
association of many stars. 

Our star formation criteria are (1) converging flow, (2) 
rapid cooling, and (3) Jeans unstable on a particle scale; 



(1) 
(2) 
(3) 



(V-w)j < 0, 

^cool < td V n • 

< is 



(6) 
(7) 
(8) 



Here i CO ol, t dyn , and t SO und are the cooling time, the dynam 
ical time of a particle, and the sound crossing time, respec- 
tively, and are expressed as 
pu 



tcool — -T ) 



td 



1 



V47rGp 
hi 

sound 5 



(9) 
(10) 

(11) 



where y, is the mean molecular weight, u is the specific ther- 
mal energy, and c s is the local sound speed. 

We assume that the star formation timescale is propor- 
tional to the dynamical timescale; 

t s f = — idyw (12) 

where the star formation parameter c is set to be 1.0. With 
c = 0.1, star formation takes place slowly, which results in 
too blue colors compared with the observation (see for 
more discussion). The star formation rate (SFR) is defined 
as 

' C -p = -cV^KGpi, (13) 



Dp, 
Dt 



where p* is the stellar density. In other words, we adopt the 
ISchmidtl (^959) law, where the SFR is the power of the gas 
fraction, and the power index is n s — 1.5. This is consistent 
with the Ha observation o f disk galaxies that implies n s — 
1.3 ± 0.3 l|Kennicuttlll989h . 

The star formation criteria are estimated with the time 
interval of At 3 {, which is set to be 2 Myr. The star formation 
probabil ity P for a g as particle forming stars during At is 
given bv iKatz) (ll992Tl as 



exp 



tdyi 



-At si 



= 1 



exp 



-cy/iwGpAtJ . (14) 



A random number between and 1 is drawn to determine 
whether the gas particle forms stars during At s f. If P is 
larger than the random number, star formation occurs. P 
is larger for higher density, and the typical density to form 
stars is higher with smaller c. For c = 1.0 and c = 0.1, most 
stars form in the region with p > 10~ 24 and 10 -22 [g cm -3 ], 
respectively. Practically, the lower limit of the gas density is 
given by c. 

When the above criteria (Eq. [6-7] and Eq.[14]) are sat- 
isfied, a part of material of the gas particle turns to a star 
particle, whic h newly forms near the gas p article. We follow 
the scheme in lNakasato fc Nomotol l)2003|) . where the initial 
mass of the star particle, mj, is derived from the integration 
of the SFR over the time interval At s {', 



p-rvhi 





c . 




^1 — exp 


-- — At s£ 

£dyn 


) 







(15) 



2.2.3 Feedback 

The evolved stars eject surrounding materials and heavy el- 
ements via stellar winds and supernova explosions. Those 
heat up, accelerate, and enrich the circumstellar and in- 
terstellar medium. High energy explosions like supernovae 
produce high temperature and low density regions in the in- 
terstellar medium. In SPH methods, the numerical accuracy 
for high density regions is much better than in mesh based 
methods, but the accuracy for low densit y regions is poorer. 
In an SPH simulation such as those of iNavarro fc White) 
( 1993), the numerical resolution (100-1000 pc) is larger than 
the typical size of supernova remnants (several tens of pc). 
Thus, because of the nature of the SPH method and the 
lack of resolution in current computing resources, it is nec- 
essary to simplify the release of the energy, momentum, and 
mass from stars. We therefore assume that energy and heavy 
elements that are ejected from a star particle are equally dis- 
tributed to the surrounding gas particles within a sphere of 
feedback radius n, which is a parameter that controls the 
mixing of heavy elements. We set rf = 1 kpc, which gives a 
good fit to the ch emical evolution in the Milky Way Galaxy 
l|Kobavashill2002h . 

In this paper, we distribute the feedb ack energy in 
purely thermal form. lNavarro fc White! lll993l) proposed that 
the energy produced by a supernova explosion is distributed 
to neighbor gas particles of the star particle mostly as a 
thermal energy and the rest is distributed as a velocity per- 
turbation to the gas particles; the fraction of energy in ki- 
netic form is given by a free parameter /kin- With /kin > 0, 
star formation efficiency is smaller, surface brightness de- 
creases at the center, and metal-rich gas blow out. If we 
adopt /kin = 0.1, the effective radius is too large and the 
metallicity gradient is too shallow (see Fig |14l in |IJ. 

The energy ejection rate E e from a star particle as a 
function of age t is 

E e (t) = ml (e e ,sw^sw(t) + e e ,iiR,u(t) + e e ,i a %a(t)) . (16) 

Ttsw, 72.li, and 72i a are the rates of SWs, SNe II, and SNe la, 
respectively, of which formulations are described in 82.2.41 

The energy of SWs from solar metallicity stars is esti- 
mated to be typi cally 0.2 x 10 erg from the observation of 
OB associations l)Abbottll982fl . Although there should be a 
mass dependence where massive stars eject larger energies, 
we adopted the typical value because of lack of observation. 
The chemical abun dance of the star significantly a ffects SWs 
as M oc Z o s ijLeitherer. Robert fc Drissenlll992ri . thus we 
include the metallicity effect for very massive stars; 



e e ,sw 



0.2 x 10 
0.2 x 10 51 



){-k) (m 2 , u <m<m u ) [erg]j(17) 
{m 2 ,i < m < ma.u) 



where Z is the metallicity of the star particle (see H2.2.4l for 
the upper and lower mass limits). 

Since there are supernovae with 10 times larger energy 
than typical supernovae, hypernovae, the energy of all su- 
pernovae should not be the same. However, since the energy 
distribution function of supernovae has not been established, 
we adopt typical values of 



e e ,n = 1.4 x 10 j1 (m,2,e < m ^ m 2. 
and 



[erg] 



(18) 
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e e ,ia = 1-3 x 10 51 {m ld>l < m ^ mi djU ) [erg] (19) 

for SN e II felinnikov et all2000l) and SNe la jNomoto et alJ 
Il984f) . respectively (see H2. 2. 41 for the upper and lower mass 
limits) . 

2.2-4 Chemical Enrichment 

A star particle is not a single star but an association of many 
stars. We assume that a star particle is in fact a simple 
stellar population, which is defined as a single generation of 
coeval and chemically homogeneous stars of various masses, 
i.e., it consists of a number of stars with various masses 
but the same age and metallicity. The mass, mass of heavy 
elements, and the spectral energy distribution of the star 
particle evolve as massive stars die. From dying stars, gas is 
ejected into interstellar medium by SWs, SNe II, and SNe 
la. 

The mass of the stars associated with each star particle 
is distributed according to an initial mass function (IMF). 
The IMF is assumed to be invariant to time and metallicity 



(m) oc m 



(20) 



which is normalized to unity at mj ^ m ^ m u . The- 
oretical arguments indicate that the IMF originates from 
fragmentation of a gas cloud almost independently of l ocal 
physics in the gas (lLow fc Lvnden-Belllll97^ : ISiUdll977ft . In 
the solar neigh borhood, the Salpeter IMF with x = 1.35 
JSalpetedll955l) is a good approximation to the star counts, 
and gives a goo d fit to many properties of disk galaxies 
jKennicutdll983h . For ellipticals, the flatter IMF is favore d 
to explain the red colors of giant ellipticals llKodamalll997l) . 
We then adopt x = 1.10 with the upper and lower masses 
of m, = 0.05M Q and m u = 120M©. 

The ejection rates of the mass and heavy element i (E m 
and E Zi ) from the star particle are expressed as 



E m (t) =m°( 



and 



E H (t) = m° (e J() sw(t) + e Zl ,n(t) + e* 4l i a (t)) . 



(21) 



(22) 



The ejection rates per mass are given by the following equa- 
tions; 



em,sw 



(1 — w m ) <t>{m) dm, 



max[m2 u , mt\ 



m 2,n 



1,,,. — I (1 — vj m ) 4>{ m ) dm, 
and 



2m, la = m C Q 7il a (t). 



(23) 



(24) 



(25) 



Time dependence is in the lower mass limit for integrals, 
the turn-off mass m t , which is the mass of the star with the 
main-sequence lifetime r m = t. For simplicity, the lifetime is 
determined as iDavid. Forman fc Joneslll990l) 



log r m = 10.0 + (-3.42 + 0.88 log m) log m, 



(26) 



which gives a little (1.5 times) longer lifetime for m > 20M© 
compared with the metallicity-dependent lifetime of Ko- 
dama & Arimoto (1997). In SWs and SNe II, stars eject the 



envelope materials outside the remnants. w m is the remnant 
mass fraction, which is the mass fraction of a black hole, a 
neutron star, or a white dwarf, depending on the initial mass 
m. For SNe la, all of the evolved He core (i.e., C+O white 
dwarfs) is ejected, and the mass of the white dwarf at the 
SN la explosion is mco = 1.38Mq. 

Heavy elements are also ejected at the following rates; 



e Zi ,sw 



(1— W m 

max[m2 iU , m±\ 



Zi <j>{m) dm, (27) 



and 



Pz,m,u 4>(m) dm (28) 

' max[m2 £ , mt] 

+ / (1 — w m — p Zi m,ll) Zi 4>{m) dm, (29) 

J max[m2 1 , m±] 



e 2 .,Ia = Pz^m, la % a (£). 



(30) 



In SWs and SNe II, the heavy elements in the envelope 
are ejected. In SNe II and SNe la, the explosive nucleosyn- 
thesis takes place, and heavy elements are newly produced. 
p Zi m,n and p Zi m,ia are the stellar yields, which are the mass 
fractions of newly produced and ejected heavy elements i, 
and are given from the supernovae nucleosynthesis model 
(Nomoto et al. 1997ab) with p Zjm ,n = for m < 10M©. 
The upper and lower limits of SNe II are 7712, u = 50M© and 
m% f i = 8Mq, respectively. For stars with 50 — 120Mq, all 
mass of He core turns to a black hole. Although some metals 
such as carbon are produced in SWs, we here neglect them 
because their contribution is much smaller than that of su- 
pernovae. The dependence of w m , p Zi m,n and p Zi m,ia on the 
stellar metallicity has not been included. 

The rates of SWs and SNe II, 7?.sw and TZn, are ob- 
tained as 



Tisw = 



and 



— 6(m) dm, 



Tin 



max [m 2,u> m t] 



max mo p , m+ 



(m) dm. 



(31) 



(32) 



For the SN la rate 7?.i a , there are several alternative 
role debate. Here we adopt an SN la model based on 
the single degene r ate sc e nario including metallicity e ffects 
llKobavashi et al] Il998l : iKobavashi. Tsuiimoto fc Nomotol 
l200dk 



m ld,u 

r 1 r. 



(m) dm 



■ <f>d(m) dm. 



(33) 
(34) 



The primary star is a C+O white dwarf (WD) formed 
from a stars with initial mass between mi Pi < = 3Mq and 
mi PtU = 8Mq — m,2,e. This SN la scenario has two types 
of secondary stars, main-sequence (MS) and red-giant (RG) 
stars. We calculate the SN la rate for each binary system 
(i.e., the MS+WD and the RG+WD system) with respec- 
tive b, mia,i, and mid tU , and combine them. The mass ranges 
of the companion stars are given by the simulation of binary 
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evolution, and are 0.9— 1.5-A/q for the RG+WD system and 
1.8 - 2.6M Q for the MS+WD system. Thus the lifetime of 
SNe la is ~ 2 — 20 Gyr and 0.5 — 1.5 Gyr, respectively. 
b is the binary parameters, which is the fraction of white 
dwarfs that eventually produce SNe la, and we adopt the 
same value which are determined from the chemical evo- 
lution in the solar neighborhood; [&rg , 6ms] = [0.02, 0.05] 
jKobavashi. Tsuiimoto fc Nomotoll2000l) . 

The photometric evolution of a star particle is identical 
to the evolution of the simple stellar population, of which 
spectra fx are taken from Kodama & Arimoto (1997) as a 
function of age t and metallicity Z. The absolute magnitude 
is calculated by using passbands to the the photometric sys- 
tem, which are the same as Kodama (1997). 

2.3 Initial Condition 

A cosmological initial condition is generated with the fol- 
lowing three steps. Throughout the paper, we set the cos- 
mological parameters of Ho = 50 km s _1 Mpc -1 , Q, m — 1.0, 
SIa = 0, and erg = 1.0. (In the following, lengths and masses 
are for Hq = 50.) First, we generate a periodic boundary 
condition with the lattice size 5 Mpc having a top-hat per- 
turbation of amplitude la or 3a in radius 1.25 Mpc by us- 
ing the COSMICS package (lBertschingerlll995l) . This pack- 
age uses the standard Zel'do vich approximation jZel'dovicbl 
Il97d : lEfstathiou et alJll985h to compute the displacements 
and velocities of the dark matter particles from a Gaussian 
random density field. The power spectrum of this density 
field represents the CDM spectrum for the above cosmologi- 
cal parameters. To gener ate constrained den sity field, a path 
integral method is used feertschinge dll987l) . We set the pa- 
rameters in COSMICS to get the density field realization 
with the desired number, mass, and the starting redshift of 
2 ~ 25. By changing the seed for the random number gener- 
ator, we obtain a different density field. Giving the seed at 
random, we obtain a random sample of initial conditions. 

Second, we pick up particles in a spherical region with 
the comoving radius of ~ 1.5 Mpc, the mass of ~ 1O 12 M0 
(baryon fraction of 0.1), and AT to t particles (the half for gas 
and the rest for dark matter). In this paper, we set two 
different resolutions; iVtot ~ 10000 and 60000. The mass of 
a dark matter particle is ~ 1.8 x 10 8 M Q and ~ 3.0 x 10 7 M Q , 
and the mass of a gas particle is ~ 2.0 x 1O 7 M0 and ~ 
3.3 x 1O 6 M , respectively. 

Third, we give the initial angular momentum to the 
system in rigid rotation because the simulated field is not 
enough large to generate tidal torque. The typical spin pa- 
rameter of a virialized halo in CDM cosmology ranges from 
0.01 to 0.1 and the median is 0.05 accor ding to the numeri- 
cal simulations fe.g. JWarren et alll992T) . We adopt the con- 
stant spin parameter A as small as ~ 0.02. With larger A, a 
spiral galaxy form if the galaxy does not undergo a major 
merger. If a major merger occurs, an elliptical galaxy form, 
even if the spin parameter is initially set to be as large as 
A ~ 0.1, Thus, we should note that such ellipticals are not 
included in our sample (see [2] for more discussion). We also 
add the corresponding Hubble velocity to the velocity field 
of the sphere since the equation of motion is integrated not 
in comoving units but in physical units. 

This scheme is similar to that in iKatzl dl992h. and 
has been used in many works (e.g., ISteinmetz fc Mulled 



Table 1 . Number of simulated galaxies 







run 


ellipticals 


dwarfs 


high resolution (TV < 


- 60000) 


13 


18 





low resolution (N ^ 


- 10000) 


59 


60 


46 


total 




72 


78 


16 



11994 iNakasato fc Nomotol Eooi iKawata fc Gibson! I2003D . 

Simulations from true CDM initial conditions excluding 
the artificial bou n dary effect were firstly carried out by 
iNavarro fc White! (j.994). We should keep in mind the 
boundary effect of our initial condition; materials in the 
outside of the simulated sphere is neglected, which makes 
the artificial cut off of the mass accretion (see for more 
discussion) . 

3 RESULTS 

3.1 EVOLUTION HISTORIES 

We simulate the chemodynamical evolution of 72 fields with 
different cosmological initial conditions (see H2.31 . It takes 
~ 10 days to simulate one field with high-resolution (the 
number of particles N ~ 60000), and 1 day with low- 
resolution (N ~ 10000). By the present time (i.e., t = 13.2 
Gyr, z = 0), in 42 cases one galaxy forms in the center of 
the field, in 16 cases one galaxy and several subgalaxies, and 
in 14 cases a few galaxies with comparable masses. We se- 
lect galaxies having stellar masses in a 20 kpc sphere larger 
than 4.5 x 1O 7 M0. Although many less-massive subgalaxies 
form, we discard them because our resolution is not enough 
to study them in detail. We summarize the number of runs 
and the resulting galaxies in Table 1. 

Figure shows the V- luminosity distributions in ±10 
kpc of 124 simulated galaxies. Luminosities are projected 
on the X-Y plane (the rotational axis is Z), and smoothed 
over 0.3 kpc. The number on the top of each image is the 
ID number in our simulation. The galaxies with the same 
number are in the same simulated field. The character "H" 
and "L" respectively denotes the high- and low- resolution, 
and the following "1" and "3" respectively denotes the la 
and 3a over-dense regions. The 18 galaxies in the left side 
are obtained with high-resolution simulations, which clearly 
show triaxial distributions. The stellar masses span in the 
range ~ 1O 8_11 M0, and cD galaxies are not included in our 
sample. Among the 124 galaxies, 78 bright galaxies are el- 
liptical galaxies, which have de Vaucouleurs' surface bright- 
ness profiles. Since the initial angular momentum is set to 
be small (spin parameter A ~ 0.02), no spirals form. The re- 
maining 46 faint galaxies are either dwarf ellipticals or dwarf 
irregular galaxies, all of which have small stellar masses as 
M < 1O 9 M . Some dwarf galaxies show bright central cores, 
while others are diffusely distributing stellar systems. 

Different galaxies undergo different evolution histories. 
The difference is seeded in the initial condition. Galaxies 
form through the successive merging of subgalaxies with 
various masses, which varies between a major merger at 
one extreme and a monolithic collapse of slowly-rotating gas 
cloud at the other. We show examples of the two cases from 
the high-resolution simulations in Figures Eland 01 which 
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Figure 1. V-luminosity distributions in ±10 kpc of 124 simulated galaxies with the stellar masses of ■ 
top of each image is ID number in the simulation. 



10 s 



Mq. The number on the 



show the evolution in ±100 kpc on the X-Z plane of dark 
matter (first lines), gas (second lines), stars (third lines), V- 
luminosity (forth lines), luminosity weighted stellar metal- 
licity (fifth lines) of galaxies that form through monolithic 
collapse and through major merger, respectively. 

Figure [5] shows the evolution of a galaxy (ID 
H3/782389913) that forms monolithically. At the beginning, 
the system expands according to the Hubble flow. The CDM 
initial fluctuation produces the structures of nodes and fila- 
ments. Gas cores form in the nodes, and stars form in the gas 
cores. The surface brightness is as high as 14 mag arcsec -2 
in rest-frame V-band. Gas rich subgalaxies merge with one 
another, and the protogalaxy coalesces at z > 3. The accre- 
tion of small subgalaxies continues till z ~ 2, and after this 
no significant event happens. There is only small amount of 
star formation at z < 0.5, and the luminosity of the galaxy 
decreases gradually toward 2 = 0. 

Figure |3] shows the evolution of a galaxy (ID 
H3/402905511) that undergoes a major merger at z ~ 2.0. 
The primary galaxy forms in a large gas core at z > 3. 
The secondary galaxy forms at the same time about 200 
kpc away, which comes to the center because of gravity, 
and merges with the primary galaxy. The secondary initially 
passes through the primary (z = 1.9), and then oscillates de- 
creasing in mass until the system is relaxed at z — 0.9. We 
should note that the edge of the dark matter distribution 



can be seen to fall onto the central object at z ~ 6.8 (Fig|5J 
and 5.0 (Fig|3J. This means that the mass accretion after 
then may be underestimated. 

We classify elliptical galaxies into the following 5 classes 
according to their merging histories. 

[1] Monolithic — Galaxies form through the assembly of 
many (> 10) gas-rich subgalaxies with the stellar masses 
of M ~ 10 9 A/q. Such assembly generally has finished by 
z ~ 3, at least by z ~ 2. The gas fractions (i.e., the mass 
ratio of gas to baryon / g = M S /(M + M S )) of merging galax- 
ies are larger than 0.5. The material of subgalaxies quietly 
accretes on the central galaxy. It is difficult to discriminate 
these subgalaxies, and this assembly looks like a monolithic 
collapse. 

[2] Assembly — Galaxies form through the assembly of 
subgalaxies with M ~ 1O 1O M0. The gas fractions / g of sub- 
galaxies are as large as 0.4. Each subgalaxy has an evolved 
core, which violently merges with the others. While the sub- 
galaxy passes through the central galaxy many times, many 
stars of the galaxy are stripped and some of them accrete 
on again. 

[3] Minor merger — The formation of the main com- 
ponent of the present-day galaxy is the same as above two 
classes, but these galaxies undergo minor merger events at 
z < 3. We define the minor merger when the stellar mass 
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Figure 2. The evolutions in ±100 kpc of dark matter (first lines), gas (second lines), stars (third lines), V-luminosity (forth lines), and 
stellar metallicity (fifth lines) of the galaxy that forms monolithically. The metallicity range is log Z/Zq = —1 to 0.4. 
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Figure 3. The same as Figure[2] but for the galaxy that undergoes a major merger at z ~ 2.0. 
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Figure 4. The star formation rates (log SFR [Mq yr -1 ]) as functions of time f [Gyr] in a present-day galaxies (r 20 kpc and \Z\ sC 100 
kpc). 



ratio of the merging galaxies / = M2/M1 (Mi ^ M2) ranges 
from ~ 0.01 to ~ 0.2. With such minor merger events, the 
surface brightness profile and the metallicity gradients are 
not affected so much. 

[4] Major merger — Galaxies undergo the major merger 
with / > 0.2 at z < 3. The mass of the primary galaxy is 
no more than five times larger than the mass of the sec- 
ondary galaxy. The redshift z ~ 3 generally corresponds 
to the galaxy formation epoch, and the main component of 
the present-day galaxy forms at z > 3. Thus, the major 
merger occurs after the most stars in the present-day galax- 
ies formed. The merger event destroys the metallicity gradi- 
ent that has existed in pre-merger galaxy in a way depending 
on the mass ratio / and the gas mass of the secondary galaxy 
(see EHl . 

[5] Multiple major mergers — Galaxies undergo a ma- 
jor merger (/ > 0.2) and one or two other mergers with 

Dwarf galaxies are classified into the following 4 classes ac- 
cording to their star formation histories. Observationally, 
the first class of galaxies are dwarf ellipticals, the others are 
dwarf irregulars. 

[Dl] Initial star burst — Galaxies form with the initial 
star burst at z > 1. In some galaxies, many supernova ex- 
plosions occur and cause the galactic winds. In other galax- 



ies, the gas is not ejected completely, but the gas density is 
so small that only few stars form at lower redshifts. Thus 
the colors are red, and these dwarfs follow the same color- 
magnitude relation as giant ellipticals. 

[D2] Continuous star formation — Galaxies grow 
though continuous star formation. After the initial star 
burst, there are the accretion of gas clumps and/or the in- 
teraction with other galaxies, which make the star formation 
continue to lower redshifts. 

[D3] Continuous star formation with recent star 
burst — The same as [D2], but a star burst occurs at the 
recent 2 — 3 Gyr. Thus the galaxy colors are blue. 

[D4] Recent star burst — Galaxies form through recent 
star bursts at z ~ 0.7. Such star bursts are induced by gas 
accretion and/or galaxy interactions. The colors are blue, 
and it is impossible to distinguish [D3] and [D4] with colors 
alone. 

The number of galaxies in each class is [1] 5 (4.0%), [2] 18 
(15%), [3] 19 (15%), [4] 25 (20%), [5] 11 (8.9%), [Dl] 20 
(16%), [D2] 13 (10%), [D3] 9 (7.3%), and [D4] 4 (3.2%). 
The percentages of non-major merger ([l]-[3]) and major 
merger galaxies ([4]-[5]) are 34% and 29%, respectively. We 
should note that the number of merger galaxies tends to be 
underestimated because of the boundary effect of our initial 
condition. 
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H3/1S9341044 HI/1342944 L3/1010791 1 L3/56971700 L3/935635703 L3b/566951 465 LI/130056237 LI /1 981 1 221 7d LI /41 799B955d LI/508323408 LI /579622B5d LI/971817698 Lib/555996139 LI b/75950 1 503h- 





H3/211903249 H1/1342944b L3/179677819 L3/638728146 L3/935635703b L3b/584285S71 L1 /1 30056237b L1/223425799 Ll/429900808 L1/50902031 3 L1/761727802 L1 b/1 5695761 9 L1 b/5859961 39b LI b/7 59301 563c 



^w: ^ 




H3/323857814 H1/1342944c L3/182721593 L3/684326117 L3/944752816 L3b/70828455 1 L1 /1 30056237c LI /223425799b LI/430769963 LI/50902031 3b LI /761 727802b LI b/1 5695761 9b L1 b/5859961 39c LI b/75930 1 563d 







H3/362982317 H1/1342944d L3/384389200 L3/69 1676562 L3/94475281 6b L3b/70828455 1 b L1/1 941 4523 L1/3B9846712 LI /430769963b L1/541 66480 L1 /761 727802c L1 b/1 5695761 9c L1b/62540932 L1 b/9 1 23B871 7 




H3/402905511 H1/3031B7717 L3/478898490 L3/706820511 L3/950724159 L3b/81 53 1 5296 L1/19414523b L1 /38984671 2b L1/47609478 LI/561948824 L1 /761 727802d L1 b/29827433S LI b/62540932b L1 b/91 238871 7b 




H3/53054973 HI/390367807 L3/491686691 L3/74644244 LZ/9"-Z-74^-, 1.31 I. Ll/19414523c LI /3S984671 2c L1/47609478b LI /561 948824b LI/785174866 LI b/29B274338b L1 b/ 6 2 5409 32c Lib/988792545 




H3/621981817 H1/985231564 L3/491 68669 1 b L3/787724242 



L3b/815315296c L1/198112217 L1/417998955 LI/485194732 L1/57962285 L1 /7851 74866b L1 b/29S274338c L1 b/749073524 L1 b/988792545b 




H3/7B2389913 HI /9B5231 564b L3/528962188 L3/7B7724242b L3b/397005886 L3b/931 089506 LI /1 981 1 221 7b LI /41 7998955b LI /4851 94732b L1/57962285b LI/789498055 Lib/343935102 LI b/749073524b 




~B 





H3/991766558 HI /985231 554c L3/563B69921 L3/6016B1B5 L3b/399652532 L3b/969939401 LI /1 981 1 221 7c LI /41 7998955c LI /4B51 94732c L1/57962285c L1/B54304394 L1 b/50540901 7 Lib/759301563 
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Figure 5. The metallicity gradients, where [M/H] is plotted against log r/r c . The thick and thin lines show the oxygen and iron 
abundance gradient, respectively. The fitting lines are shown with the dotted lines. The dotted vertical lines show the maximum radii 
used in the fitting. 



Figure 0| shows the star formation rates (log SFR [Mq 
yr -1 ]) as a function of time t [Gyr], which are derived from 
the ages and masses of stars that belong to the galaxy at 
present. Thus, these rates are not for stars formed in the 
region that is defined as the galaxy at each redshift, but 
for stars that formed anywhere and are today part of the 
galaxy. This definition of the SFR will be the same as in 
any observation that estimates the star formation history of 
a galaxy from its stellar populations. As clearly shown, all 
elliptical galaxies form with an initial star burst at z > 2, 
whereas dwarf galaxies undergo relatively continuous star 
formation. The SFR decreases because the gas is exhausted 
in the galaxy. The secondary star burst is induced by the ac- 
cretion of gas clumps and/or the merging of gas-rich galax- 
ies. Not all merging events induce a secondary star burst; 
the fraction of merging events that induce such a star burst 
is about 10%. For all ellipticals, the initial star burst is al- 
ways larger than the secondary one. The typical timescale 
of initial star burst is found to be 1 — 2 Gyr, which is much 
longer than the 0.1 Gyr that is commonly ado pted in the 
one-zone model (e.g.. lKodama fc Arimotolll99'it) . Such SFR 
is due to the artificial cut-off of mass accretion caused by 
our initial condition, but is required from the observation of 
ellipticals as summarized in 

We should note that star formation has not completely 



stopped at present in most elliptical galaxies. Galactic winds 
are hard to generate in our simulated giant ellipticals (see 
for more discussion) . At the center of the present-day 
galaxy, the dynamical potential is so deep that the gas den- 
sity is high. In such regions, super metal-rich stars (Z ~ 
IQZq) form in the simulation. In several dwarf galaxies with 
M < 1O 9 M (e.g., ID Ll/223425799b and Ll/430769963), 
weak galactic winds can be seen even if /ki n = 0. However, 
enough gas is heated up and blows away gradually by the 
input of thermal energy of supernovae. The global gas frac- 
tion of the simulated fields spans over 50 — 90%, and the gas 
fraction in a galaxy (a sphere of 2r c ) is 1 — 10% and 30 — 80% 
for giants and dwarfs, respectively. The fraction of heavy el- 
ements locked into stars in a galaxy sp ans 20 — 40% which 
is consistent with the observation (e.g.. iRenz ini 2002). 

3.2 METALLICITY GRADIENTS 

Galaxies are observed in projection on the sky. In order to 
compare the simulated results with observation, foreground 
and background particles should be excluded properly. We 
define a galaxy as the projection of \Z\ ^ 100 kpc on the X- 
Y plane. Below radius r means a projected radius. We then 
derive the effective radius r c by fitting a de Vaucouleurs' law 
to the surface brightness profile. For the simulated galaxies, 
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Figure 6. Mctallicity gradients versus (a,d) stellar masses in r e , (b,e) total luminosities, and (c,f) mean stellar metallicities in r c . The 
upper (a-c) and lower panels (d-f) show the gradients of oxygen and iron abundances, respectively. All points are for the simulated 
galaxies, and the large points denote galaxies simulated with high-resolution. The shape of symbols show the merging histories for 
elliptical galaxies and star formation histories for dwarf galaxies; [1] monolithic (red filled circles), [2] assembly (yellow filled squares), 
[3] minor merger (green filled triangles), [4] major merger (cyan open squares), and [5] multiple major merger (blue open circles); [Dl] 
initial star burst (asterisks), [D2] continuous star formation (crosses), [D3] continuous star formation with recent star burst (plus), and 
[D4] recent star burst (three-pointed stars). 



we exclude the central region with r — — 1 kpc and — 2 
kpc when fitting simulations with high- and low- resolution, 
respectively, because the surface brightness is smeared out 
due to the softening of the gravity. In the outer region, to 
derive proper effective radii in the simulation, we use large 
enough regions such as 40 — 80 kpc. In some case that galax- 
ies show a local excess in the surface brightness because of 
the existence of the satellite subgalaxies, we only use the 
inner part. 

Observationally, metallicity gradients are derived from 
the line-strength gradients, and line-strength is converted to 



metallicity using an index-meta llicity relation derived from 
spectral synthesis models (e.g.. iKodama fc Arimotolll99^ . 
Thus the observed metallicity is the luminosity weighted 
metallicity. In the observational data, the gradients are 
smeared out due to poor seeing conditions at galaxy cen- 
ters with log r/r c < — 1.5, and in the outer regions with 
> 2r e errors arising from the sky subtraction give poor fits. 
Kobayashi & Arimoto (1999) excluded these regions from 
the fitting. In the simulation, we provide the metallicity 
weighted by V-luminosity, because the absorption indices 
usually observed such as Mgg and Fei are in the V-band. 
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Since we find no significant difference in the metallicity gra- 
dients at r = 1 — 2 kpc for the high- and low-resolutions, 
the innermost boundary of the fitting region is set at 1 kpc 
for both resolutions. This value is larger than the radius of 
log r/r c ~ —1.5. The outer boundary is set at 2r e for most 
galaxies, and at 10 kpc for galaxies in the case with small 
r c . The contribution of satellite galaxies is also excluded. 

Figure [S] shows the metallicity gradients of 124 galax- 
ies, where [M/H] is plotted against log r/r c . The thick and 
thin lines show the oxygen and iron abundance gradient, re- 
spectively. The fitting lines are shown with the dotted lines. 
The dotted vertical lines show the maximum radii used in 
the fitting. As is clearly seen, the metallicity gradients are 
various, some are steep and others flat. 

3.3 RELATION AGAINST MASS? 

Figure HJ shows the V-luminosity weighted metallicity gra- 
dients versus stellar mass within r e (a and d), total lumi- 
nosities derived from the de Vaucouleurs' law (b and e), and 
luminosity weighted mean stellar metallicities in r c (c and 
f). The upper (a-c) and lower panels (d-f) show the gradi- 
ents of oxygen and iron abundances, respectively. All points 
are for the simulated galaxies, and the large points denote 
galaxies simulated with high-resolution. The symbols show 
the merging histories for elliptical galaxies and star forma- 
tion histories for dwarf galaxies; [1] monolithic (filled cir- 
cles), [2] assembly (filled squares), [3] minor merger (filled 
triangles), [4] major merger (open squares), and [5] multiple 
major merger (open circles); [Dl] initial star burst (aster- 
isks), [D2] continuous star formation (crosses), [D3] contin- 
uous star formation with recent star burst (plus), and [D4] 
recent star burst (three-pointed stars). 

The remarkable result is that there is no correlation be- 
tween the gradients and masses or luminosities in Figure [5] 
The lack of these relations has already been noticed in the 
observat ional data, and the origin o f the scatter has been 
argued llKobavashi fc Arimotol I1999T) . It seems that more 
metal-rich simulated ellipticals tend to have steeper gradi- 
ents, but it is hard to find any relation if simulated dwarfs 
are included. The origin of scatter is clearly shown with the 
symbols; the galaxies that form monolithically (filled circles 
and squares) have steeper gradients, and the galaxies that 
undergo major mergers (open circles and squares) have shal- 
lower gradients. Therefore, we conclude that the metallicity 
gradients do not depend on the galaxy mass, and the vari- 
ety of the gradients stems from the difference in the merging 
history. 

The distributions of metallicity gradients for [1] mono- 
lithic collapse, [2] assembly, and [3] minor merger are very 
similar. The minor mergers do not affect so much the metal- 
licity gradients. The distributions for [4] major merger and 
[5] multiple major mergers are also similar. Then we joint 
the 5 classes into 2 large classes: [A] non-major merger 
galaxy including [l]-[3], and [B] major merger galaxy includ- 
ing [4] and [5]. Figure |7| show the histograms of the metal- 
licity gradients for the 2 classes using metallicity Z (top 
panel), oxygen (middle panel), and iron abundance (bottom 
panel). The distributions for [A] non-major merger (gray 
area) and [B] major merger (hatched area) galaxies are dif- 
ferent. The typical gradients for non-major merger and ma- 
jor merger galaxies are A log Z/ A log r ~ —0.30 and —0.22, 




-0.8 -0.6 -0.4 -0.2 

AlogZ/Alogr 




-0.8 -0.6 -0.4 -0.2 

Alog[0/H]/Alogr 




-0.8 -0.6 -0.4 -0.2 



Alog[Fe/H]/Alogr 

Figure 7. The histograms of the gradients for metallicity Z 
(upper panel), oxygen (middle panel), and iron abundance (lower 
panel). The gray and hatched area are for non-major merger (red) 
and major merger (blue) galaxies, respectively. The thick lines 
show the observation with Mg2 (upper and middle panels) or Fei 
(lower panel) index from Kobayashi & Arimoto (1999). 

A[0/H]/Alogr ~ -0.30 and -0.18, and A[Fe/H]/A log r ~ 
—0.45 and —0.38, respectively. The galaxies with gradi- 
ents as steep as AlogZ/Alogr ^ —0.35 are all non-major 
merger galaxies. 

The thick lines show the observed metallicity gradient 
distributions for Mg2 (top and middle panels) and Fei (bot- 
tom panel) using data from Kobayashi & Arimoto (1999). 
(The reason to compare oxygen gradients with the Mg2 ob- 
servation is that the yields relative to solar value are almost 
the same between O and Mg.) The simulated Z gradients 
are consistent with the Mg2 observation both with the mean 
value of A log Z/A log r ~ —0.3 and the dispersion of ±0.2. 
The oxygen gradients are shallower by 0.05 dex than the Mg2 
observation. This is because the Mg2 index is not a good in- 
dic ator for magnesium abundance, but rather for met allic- 
ity dTripicco fc BelJll995t iKobavashi fc Arimotolll99^ . For 
iron, the simulated Fe gradients are steeper than the ob- 
servation, even if one takes into account the much larger 
observational errors of Fei than Mg2. This is because in the 
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simulation the star formation does not terminate completely 
in the galaxy center (see jgj. 

From the observed metallicity gradients, we can es- 
timate the fraction non-major merger and major merger 
galaxies. Using the Mg2 observation, we derive the num- 
bers of nearby elliptical galaxies with steep gradients of 
A log Z/ A log r ^ —0.3 and flat gradients of > —0.3 as 25 
and 21, respectively. If we take alternative threshold value 
of A log Zj A log r = —0.25, the numbers become 30 and 16, 
respectively. Therefore, the fraction of non-major merger 
galaxies is half or two third. Even if we take account of 
some problems involved in our initial condition (see 0for 
the detail), there exist non-major merger galaxies and ma- 
jor merger galaxies half and half. The observed variation 
in the metallicity gradients cannot be explained by either 
the monolithic collapse only or the major merger only. It is 
well reproduced in the present model where both formation 
processes arise under the CDM scheme. 




T - 



t [Gyr] 

Figure 8. The evolution of metallicity gradients for the mono- 
lithic collapse case. The thick and thin lines are for the oxygen and 
iron abundances, respectively. The solid and dotted lines are the 
gradients weighted by V-luminosities and masses, respectively. 



3.4 EVOLUTION OF GRADIENTS VIA 
MERGERS 

The spatial distribution of metallicity is already shown in 
the last lines of Figures ([1] the monolithic case) and |H] 
([4] the major merger case). The metallicity is enhanced in 
the high density region, and radial metallicity gradients ap- 
pear in protogalaxies at z ~ 5. After the major merger at 
z ~ 2.0, many metal-rich stars move to the outer region of 
the galaxy. In the monolithic case, metal-rich region is con- 
centrated, while metal-rich region extends over ±10 kpc in 
the major merger case. The present metallicity gradients are 
A log Z/Zq I A log r = -0.46 and -0.19, respectively. 

Figures |S] and |U] show the time evolutions of metallicity 
gradients for [1] the monolithic case and [5] the multiple ma- 
jor merger case, respectively. The thick and thin lines are for 
the oxygen and iron abundances, respectively. The solid and 
dotted lines indicate the gradients weighted by V-luminosity 
and mass, respectively. The gradients are underestimated if 
the metallicity is weighted by mass. 

Figure |H| shows the gradient evolution of ID 
H3/782389913. This galaxy forms through the assembly of 
3 galaxies with the stellar masses of 8x 10 9 M Q , 4x 10 9 M Q , 
and 4 x 1O 9 M0, and thus the initial gradient is not so steep, 
A[Fc/H]/A log r ~ —0.9. This assembly has finished by 
z ~ 3.7 (t ~ 1.2 Gyr). Then the galaxy evolves quietly, and 
the metallicity gradient is kept nearly constant. At z ~ 0.1 
(t ~ 10 Gyr), a small galaxy with ~ 1O 8 M0 accretes. The 
small star formation is induced, which changes the gradient 
by 0.1 dex temporarily. 

Figure [5] shows the gradient evolution of ID 
H3/390367807. This galaxy undergoes the major mergers 
three times at z ~ 1.3,0.9, and 0.5 (t ~ 3.7,5.1, and 6.9 
Gyr). The mass ratios are / = 0.47,0.16, and 0.24, re- 
spectively. The initial star burst produces the metallicity 
gradient of A[Fe/H]/A log r ~ —1.2. The major merger 
event is so violent that the gradient decreases by 0.7 dex 
to be A[Fe/H]/Alogr ~ -0.5. The second and third 
merger events are not so active, and the gradient becomes 
A[Fe/H]/A log r ~ —0.4 and —0.3 at the end of the merger 
event, respectively. 

The initial gradient is determined during the initial star 
burst at z >3. The gradient is steeper in the case of quiescent 
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t [Gyr] 

Figure 9. The same as Figure 03 but for the major merger case. 
Mergers takes place at t ~ 3.7, 5.1, and 6.9 Gyr. 

gas accretion, and is shallower in the case of violent assembly 
of subgalaxies. As a result, the initial gradients span from 
A[Fc/H]/A log r ~ —1.5 to —1.0. Since the gradients evolve 
a lot after this time, the initial gradients cannot be inferred 
from the present gradients. 

The metallicity gradients become definitely shallower 
when galaxies merge. However, the formation and destruc- 
tion of gradients via mergers are complicated. To find out 
some rules, we examine the physical conditions at the 151 
merging events that occur in our simulations of giant galax- 
ies. Figure tTfjl shows mass ratios / = M2/M1 (Mi ^ M2) (a), 
ratios of gas mass M gl 2/M Sl i (b), gas fractions of primary 
galaxies / gj i (c), and those of secondary galaxies f St % (d) 
against redshifts z. The averages of / and / gj i in each red- 
shift bin decrease toward z — because the stellar masses 
of primary galaxies increase. The average of M gl 2/M g ,i is 
larger than that of /, and is constant of ~ 0.5. / g ,2 is larger 
than / g ,i, and spans from to 1. These are because star 
formation strongly depends on gas density and takes place 
slowly in the less-massive secondary galaxies. At z < 1, some 
major mergers (/ ,>0.2) have M gl 2/M gl i ,>0.5, which means 
that the secondary galaxies have comparable mass and gas 
content to the primaries. These induce strong star forma- 
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Figure 10. (a) Mass ratios / = M2/M1 (Mi M2), (b) ratios of gas mass Mgn/Mg 1, and (c,d) gas fractions (/ g ,i and /g 2) of 
the primary and secondary galaxies against rcdshifts z. For the symbols, see the caption of Figure The points surrounded by large 
circles/triangles indicate that strong/moderate star formation is induced. 



tion (surrounded by circles). Relatively weaker star forma- 
tion (surrounded by triangles) is induced by mergers with 
/ g ,2 > 0.5, which means that the stellar mass of the sec- 
ondary is small and the secondary galaxy behave like a gas 
cloud. 

Figure 1111 shows the change of gradient, 
A(A[Fe/H]/A log r), caused by a merging event at z < 3 
against various physical conditions: / (a), M gl 2/jWg,i (b), 
/ gl i (c), and / gl 2 (d). The gradient change strongly depends 
on the mass ratio /. The mergers with / > 0.2 decrease the 
gradients at least by 0.5 dex. This is the reason why we 
have defined major mergers as those with / > 0.2. A major 
merger changes the orbits of stars. The stars that are in 
the center and have high metallicities are able to move to 
the outer region of the galaxy. However, if the secondary 



galaxy contains as much gas as the primary galaxy (i.e., 
Af gl 2/Mg,i >^ 0.5) and strong star formation is induced 
(surrounded by circles), then the gradient change stays 
smaller than 0.5 dex. This is because such star formation 
takes place after the gas falls into the center of the primary 
galaxy and increases the metallicity at the center, resulting 
in a small gradient change. Sometimes merging of a gas rich 
galaxy (/ g ,2 > 0.5) also induces moderate star formation 
(surrounded by triangles). In this case, star formation takes 
place in the outer region of the galaxy, and the gradient 
change can be as large as ~ 0.5 dex, even if / ~ 0. In some 
cases without a merging event, a similar star formation is 
induced by late gas accretion, and the metallicity gradient 
gradually becomes shallower. 

By showing the similar figure in Ivan Albadal lll982Tl . we 
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Figure 11. The gradient change A(A[Fe/H]/A log r) before and after the merging event at z < 3 against (a) mass ratios /, (b) ratios 
of gas mass Mgi/Mgi, and (c,d) gas fractions (/ g ,i and fg,2). For the symbols, see the caption of Figure l6l The points surrounded by 
large circles/triangles indicate that strong/moderate star formation is induced. 



argue that dynamical information on the orbits of N-body 
particles is not fully wiped out by a merger, but fairly lost to 
change the metallicity gradient. Figures E21 and IT51 show the 
energies E = u 2 /2 + $ of particles before and after a merging 
event. Large E means that a particle has an extended orbit. 
Figure Hoi shows the energies at 3.5 and 4.9 Gyr for the first 
major merger at ~ 3.7 Gyr of the galaxy ID H3/390367807 
that undergoes a triple major merger (see Figure |5] for the 
gradient evolution) . Figure 1121 shows the same but for the 
non-merger galaxy ID H3/782389913. Clearly, the dispersion 
for the merger galaxy is much larger than the non-merger 
galaxy. 

In summary, whether the merging event changes the 
metallicity gradient is mainly influenced by two factors; 1) 
the mass ratio of the merging galaxies / and 2) the induced 



star formation. The evolution of the gradients via merging 
events is based on the following three processes: i) Destruc- 
tion by mergers to an extent dependent on /. ii) Regener- 
ation due to the central star formation induced at a rate 
dependent on M gl 2/M g> i. iii) Passive evolution as star for- 
mation is induced in the outer regions at a rate dependent 
on / gi 2- We note that the gradient change is underestimated 
if the metallicity is weighted not by luminosity as the obser- 
vation, but by mass as in previous simulations. 



4 DISCUSSION 

The global properties of elliptical galaxies depend mainly on 
their masses, while their metallicity gradients are much af- 
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E(t=3.5) 

Figure 12. The energies E = v 2 /2 + $ of particles in the non- 
merger galaxy. 




E(t=3.5) 

Figure 13. The energies of particles before and after the major 
merger at ~ 3.7 Gyr. 

fected by their merging history. Merging histories can thus, 
in principle, be inferred from the observed metallicity gra- 
dients of present-day galaxies. The dispersion in metallicity 
gradients for galaxies with similar merging histories is not 
small (Fig. EJ , and so it may be difficult to determine the 
merging history of an individual galaxy from its metallic- 
ity gradient. However, it should be possible to estimate the 
fractions of non-merger and merger galaxies by using the 
statistics of metallicity gradients. For example, if the frac- 
tions are estimated for field and cluster galaxies, it should 
provide information about environmental effects on galaxy 
formation. Even allowing for the uncertainties in available 
observations of nearby galaxies, our predicted difference be- 
tween gradients of non-merger and merger galaxies is large 
enough to be detected. 

We should note that there are two problems in the sim- 
ulated galaxies, which cannot be improved by changing pa- 



rameters in the present model. One is that galactic winds do 
not occur in large galaxies in the simulations, and star for- 
mation never terminates completely. The colors of our galax- 
ies thus tend to be bluer than the observation, and the B-V 
color distribution extends to 0.4 mag bluer. In our galaxy 
centers, the color gradients tend to have the opposite slope 
to the metallicity gradients, because the stellar populations 
there are young. Such late star formation takes place in the 
central 1 kpc, and this region is carefully excluded when es- 
timating metallicity gradients in this paper. This problem 
arises from our SPH method and our feedback scheme. If we 
include the kinetic feedback with /km > 0, surface bright- 
ness decreases at the center, and metal-rich gas blow out. 
However, these results in too large effective radii and too 
shallow metallicity gradients to meet the observation. Figure 
1141 shows the surface brightness profiles and metallicity gra- 
dients with (dash-dotted line) and without (solid line) /kin 
using the same initial condition. With /ki n = 0.1, the surface 
brightness decreases by ~ 2 mag at the center, which results 
in r c = 8.1 kpc (r c = 4.3 kpc for /kin = 0). The metallicity 
gradient is almost flat within r c . To describe galactic winds 
with the kinetic feedback, star formation scheme should be 
modified simultaneously. 

We note that our predictions for oxygen abundance are 
not changed so much, whereas the iron abundance should 
decrease especially at the galaxy center. In Figure |S| there 
appear to be a relation between the iron abundance and iron 
gradients. This is because late star formation at the galaxy 
center increases the iron abundance and steepens iron gra- 
dients at the same time. Such a relation can disappear if 
the galaxy model is improved. In Figure |7| the mean value 
of simulated iron gradients is steeper than observed. This 
is also caused by the same reason. If the galaxy model is 
improved so as to blow a galactic wind, iron gradients for 
non-merger galaxies will become the similar to their oxy- 
gen gradients. However, iron gradients for merger galaxies 
will be still steeper than oxygen gradients. The difference 
of iron and oxygen gradients, i.e., [O/Fe] gradients, may be 
the best indicator of the formation history of the galaxy. To 
discuss [O/Fe], observational data of high-enough S/N is re- 
quired, and the spectral population synthesis models should 
be prepared for various [O/Fe]. 

The other problem is that simulated galaxies are more 
extended, and thus there is an offset in the radius-magnitude 
relation. This problem is not in magnitude, because the 
Faber- Jackson relation can be reproduced, and the galaxy 
with given mass has reasonable luminosity. While the stel- 
lar mass-to-light ratio is consistent with the observation as 
M/Lb ~ 5 — 8, the total mass-to-light ratio is as large as 
(A/dm + Mbaryon)/£ ~ 10 — 100 (within a sphere of 2r c ), 
which may be too large. The baryon fraction increases to 
~ 0.5 at the center (r < 1 kpc) because of the stellar con- 
centration, but is as small as ~ 0.2 in r < 2r e . Baryon does 
not fully dominate even at the galaxy center in our sim- 
ulated galaxies, and the stellar concentration seems to be 
not enough compared with the dark matter concentration. 
These might be due to cosmological parameters, dynamical 
friction, star formation parameter, and the limited field size 
of initial condition. As well kn own with the zero-point off - 
set in the Tully- Fisher relation dNavarro fc Steinmetzl2 000). 
such high dark matter concentration should arise from the 
standard CDM cosmology that is adopted in this paper, and 
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Figure 14. Surface brightness profiles (upper panel) and metal- 
licity gradients (lower panel) with kinetic feedback (dash-dotted 
line, = 0.1), with longer timescale of the star formation (dot- 
ted line, c = 0.1), and wider region of initial condition (2.5 Mpc). 



should d ecrease to some extent with t he A-CDM cosmology. 
However. iNavarro fc Steinmetzl i2000l) argued that the zero- 
point offset remains even with the A-CDM cosmology, which 
requires substantial revision to the CDM scenario such as as, 
a tilt in the primordial power spectrum, or hot dark matter. 
The dynamical friction may be effective because the mass 
of a dark matter particle is ten times larger than th e gas 
and stellar particles (but see lSteinmetz fc Whiteil99'i) . Too 
large effective radius suggests that the star formation takes 
place too early before the gas accretes towards the center. 
This can be solved by changing the star formation timescale, 
i.e., reducing the star formation parameter c, although col- 
ors become bluer. As shown in Figure tm with c = 0.1 (dot- 
ted line), the surface brightness increase by 1 mag at the 
center, which results in smaller effective radius as r c = 3.2 
kpc. The metallicity increases because of longer timescale of 
star formation, but the gradient dose not change so much. 
As mentioned in H3.2I the surface brightness profiles of our 
simulated galaxies well follow the de-Vaucouleurs law, but is 
smeared out at the center, and the central surface brightness 
(r < 1 kpc) is slightly smaller than the de-Vaucouleurs fit 
(see Fia ll4t . Although the resolution is not enough to dis- 
cuss, this central stella r concentration may be sma ller than 
other simulations (e.g.. ISpringel fc Hernauistll2003l) . can be 
as large with c = 0.1, and should be affected by the differ- 
ence in the star formation and feedback schemes. 

The initial condition is most important. The edge of the 
simulated fields actually falls into a galaxy at z ~ 2 — 4 de- 
pending on the position of the galaxy in the field. The galaxy 
mass M200 (the mass with higher density than 200p cr it) stops 



increasing at z ~ 1.5 — 3. This means that there is an ar- 
tificial cutoff of mass accretion around this redshift, which 
results in the cutoff of star formation rate in Figure 2] In 
a wider simulation with a radius of ~ 2.5 Mpc (dashed-line 
in Fig |14H instead of ~ 1.5 Mpc in the series of our simu- 
lation, the star formation continues longer, and thus colors 
tend to be too blue. The total luminosity becomes ~ 2 times 
larger, but still ~ 3 times smaller than the observation. With 
much wider simulation, the luminosity difference could be 
improved, but the color inconsistency become much larger, 
and such galaxy is no more elliptical. The star formation 
and feedback schemes need to be mod ified accordingly, such 
as so me recent works have tried (e.g.. lSpringel fc Hernauisti 
2003). In this paper, the mass accretion and star formation 
are truncated artificially by the initial condition. However, 
in observed ellipticals, star formation should be truncated 
at z ~ 2 by some process; tidal stripping, effects of active 
galactic nuclei, and so on. This epoch may correspond to 
the epoch when the galaxy falls into a cluster, because the 
galaxy moves too fast to arise mass accretion. 

If we include the contribution of the outside of the simu- 
lated field, late accretion and late merging events should in- 
crease. Both make the metallicity gradient shallow. As noted 
in H2.3I fields with larger spin parameter A are not included 
in our simulation, and an elliptical can form through major 
mergers of disk galaxies in some of such fields. However, as 
shown in Figure |7| the predicted metallicity gradients are 
already in good agreement with the observation, and the in- 
crease of merger galaxies conflicts with the observation. If 
star formation and feedback schemes improved, evolution of 
metallicity gradients may also be affected at some extent. We 
should note, however, that the following conclusions will not 
change; major merger make the metallicity gradients shal- 
low, and there should be non-merger ellipticals to explain 
the observed steep gradients. 



5 CONCLUSIONS 

We study the formation and chemodynamical evolution of 
galaxies with our GRAPE-SPH chemodynamical model that 
includes various physical processes associated with the for- 
mation of stellar systems; radiative cooling, star formation, 
feedback from SNe II, SNe la, and SWs, and chemical enrich- 
ment. We simulate 72 slowly-rotating spherical fields (spin 
parameter A ~ 0.02), and obtain 124 galaxies (78 ellipti- 
cals and 46 dwarfs) from the CDM initial fluctuation. All 
simulated galaxies have the de Vaucouleurs' surface bright- 
ness profiles, and are therefore elliptical galaxies. Most stars 
in ellipticals form during the initial star burst at z > 2, 
while dwarfs undergo relatively continuous star formation. 
In our scenario, galaxies form through the successive merg- 
ing of sub-galaxies. The merging history is various and the 
difference is seeded in the initial conditions. In some cases, 
galaxies form through the assembly of gas rich small galax- 
ies, and the process looks like a monolithic collapse. In other 
cases, the evolved galaxies undergo major merger of galax- 
ies. Major mergers are defined as those with the mass ratios 
of the primary and secondary galaxies being / ,> 0.2. 

We examine the physical conditions during 151 merging 
events that occur in our simulation. Whether the merging 
event changes the metallicity gradient is mainly influenced 
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by two factors; the mass ratio of the merging galaxies / 
and the induced star formation. The basic processes of the 
formation and evolution of the gradients are summarized 
below: 

• Formation of initial gradients — The initial gradient is 
determined from the initial star burst at z > 3. The gradient 
is steeper in the case of quiescent gas accretion, and is shal- 
lower in the case of violent assembly of subgalaxies. As a re- 
sult, the initial gradients span from A[Fe/H]/A log r — —1.5 
to -1.0. 

• Destruction by mergers — The major merger changes 
the orbits of stars. The metal-rich stars at the center are 
able to move to the outer region of the galaxy. The gradient 
change is determined mainly from the mass ratio of merging 
galaxies /. With larger /, the gradients become shallower. 
If the mass ratio of merging galaxies is larger than / ~ 0.2, 
the gradient change is larger than ~ 0.5 dex. 

• Regeneration due to the induced star formation — If 
the ratio of gas mass is as large as M &t 2/M g ,i > 0.5, strong 
star formation is induced at the center of the primary galaxy, 
and the gradient change is smaller than ~ 0.5 dex. 

• Passive evolution — If the gas fraction of the secondary 
galaxy is larger than / g ,2 ~ 0.5, moderate star formation is 
induced in the outer region of the primary galaxy, and the 
gradient change becomes as large as ~ 0.5 dex, even if / ~ 0. 
In some case without merging event, if the similar star for- 
mation is induced by the late gas accretion, the metallicity 
gradient gradually becomes shallower. 

We succeed in reproducing the observations of metallic- 
ity gradients and finding the origin of the variety of internal 
structures. From the distribution functions of the gradients 
for different merging histories, we discuss the origin of ellip- 
tical galaxies. 

• The average metallicity gradient is A log Z/ log r ~ 
—0.3 and the dispersion is ±0.2, which are both consistent 
with observations of Mg2 gradients. 

• No correlation is produced between gradients and 
masses. The metallicity gradients do not depend on the 
galaxy mass, and the variety of the gradients stems from 
the difference in the merging histories; galaxies that form 
monolithically have steeper gradients, while galaxies that 
undergo major mergers have shallower gradients. 

• The metallicity gradient distributions for [A] non-major 
merger ([l]-[3]) and [B] major merger galaxies ([4] and [5]) 
are quite different. The typical gradients for non-major 
merger and major merger galaxies are A log Z/A log r ~ 
—0.3 and —0.2, respectively. Simulated galaxies with the gra- 
dients steeper than —0.35 are all non-major merger galaxies. 

The global properties of elliptical galaxies depend 
mainly on their masses, while their metallicity gradients are 
much affected by their merging history. A major merger 
makes the gradient shallower. Therefore, merging histo- 
ries can be inferred from the observed metallicity gradients 
of present-day galaxies. Available observations for nearby 
galaxies suggest that there exist non-major merger galax- 
ies and major merger galaxies half and half. The observed 
variation in the metallicity gradients cannot be explained 
by either monolithic collapse or by major merger alone. In- 
stead, it is well reproduced in the present model in which 
both formation processes arise under the CDM scheme. 
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